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ABSTRACT 

We present tomography of the circum-galactic metal distribution at redshift 1.7 to 4.5 derived 
from echellete spectroscopy of binary quasars. We find C IV systems at similar redshifts in paired 
sightlines more often than expected for sightline-independent redshifts. As the separation of the 
sightlines increases from 36 kpc to 907 kpc, the amplitude of this clustering decreases. At the largest 
separations, the C IV systems cluster similar to Lyman-break galaxies (Adelberger etal. 2005a). The 
C IV systems are significantly less correlated than these galaxies, however, at separations less than 
Ri ~ 0.42 ± 0.15 h^ 1 comoving Mpc. Measured in real space, i.e., transverse to the sightlines, this 
length scale is significantly smaller than the break scale estimated from the line-of-sight correlation 
function in redshift space (Scannapieco etal. 2006a). Using a simple model, we interpret the new real- 
space measurement as an indication of the typical physical size of enriched regions. We adopt this size 
for enriched regions and fit the redshift-spacc distortion in the line-of-sight correlation function. The 
fitted velocity kick is consistent with the peculiar velocity of galaxies as determined by the underlying 
mass distribution and places an upper limit on the average outflow (or inflow) speed of metals. The 
implied time scale for dispersing metals is larger than the typical stellar ages of Lyman-break galaxies 
(Shapley etal. 2001), and we argue that enrichment by galaxies at z > 4.3 played a greater role in 
dispersing metals. To further constrain the growth of enriched regions, we discuss empirical constraints 
on the evolution of the C IV correlation function with cosmic time. This study demonstrates the 
potential of tomography for measuring the metal enrichment history of the circum-galactic medium. 

Subject headings: cosmology: miscellaneous - galaxies: abundances - galaxies: halos - galaxies: 
high-redshift - intergalactic medium - quasars: absorption lines 
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1. INTRODUCTION 

In essentially every variant of Big Bang nucleosynthe- 
sis, the primordial universe is composed mostly of hydro- 
gen and helium, with trace amounts of the light elements 
Li, Be, and B. Heavier elements, the metals, are produced 
later by stellar nucleosynthesis, and their mass traces 
the quantity of stars formed. Since it seems reasonable 
to expect star formation to be confined to galaxies, the 
discovery of highly ionized carbon (and later highly ion- 
ized Si and O) in the tenuous intergalactic hydrogen gas 
between galaxies came as a considerable surprise (e.g., 
Meyer & York 1987; Lu 1991; Songaila & Cowie 1996). 
Since these early observations, it has become well estab- 
lished that some enrichment persists to low Lya col- 
umn densities, suggesting at least some filamentary in- 
tergalactic structures have been enriched (Ellison et al. 
2000; Schaye et al. 2003). These metals are thought 
to be ejected from galaxies or protogalaxies by power- 
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ful galactic winds. These winds circulate metals through 
a circum-galactic medium (CGM). The formation and 
evolution of this circum-galactic medium have not been 
directly measured. 

The history of circum-galactic metal enrichment di- 
rectly impacts the metallicity of the gas accreted by all 
galaxies, the properties of the first generation of stars, 
and the galactic mass-metallicity relation. Suggestions 
for when this enrichment took place span a broad range 
from before the formation of most galaxies by explosions 
of an early generation of very massive stars (Madau, Fer- 
rara, & Rees 2001; Qian, Sargent & Wasserburg 2002), 
to during or just after reionization at z w 5 — 10 by low- 
mass galaxies (Gnedin & Ostriker 1997; Ferrara et al. 
2000; Scannapieco, Ferrara, & Madau 2002), or even sig- 
nificantly later, redshifts 2-4 (Adelberger etal. 2005b), 
when the cosmic star formation rate peaked. Observa- 
tions that constrain the period of metal dispersal are 
therefore of great interest. 

Previously, the evolution of the carbon mass den- 
sity has provided the primary empirical constraint 
on intergalactic enrichment. The carbon mass 

density inferred from the detections of intervening 
C IV AA1548.204, 1550.781, however, depends sensitively 
on the gas density and the ambient radiation field. In cos- 
mological models the fraction of carbon in the third ion- 
ization state reaches its maximum value near z sa 5 (Op- 
penheimer etal. 2009). Beyond z> 5, both the decreas- 
ing ionization fraction and the decreasing C mass den- 
sity contribute to the decline in the density of C IV sys- 
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terns (Becker etal. 2009; Ryan- Weber etal. 2009). At 
z < 5, where the density of C IV systems remains roughly 
constant (Songaila 2001), substantial IGM enrichment 
may be masked by the decline in the ionization frac- 
tion n(C +3 )/n(C) with cosmic time (Cooksey et al. 2009; 
D'Odoricoetal. 2010). 

Compared to mass density measurements, measuring 
the size of enriched regions would be less sensitive to 
the evolution of the gas density and ambient radiation 
field. Size measurements require only metal markers, like 
C IV detections, in multiple sightlines. They can directly 
constrain the period of metal dispersal because larger 
enriched regions at a given cosmic epoch require earlier 
enrichment (Scannapeico 2005). Size measurements for 
C IV systems are complicated by both the clustering of 
galaxies and the degeneracy between Doppler shifts and 
line-of-sight distance. By comparing previous strategies 
for measuring the size of enriched regions and introducing 
a new method in this paper, we aim to establish the via- 
bility of size measurements and demonstrate their utility 
for constraining the history of circum-galactic metal en- 
richment. 

The physical distance winds travel is thought to be rel- 
atively independent of galaxy mass and redshift (see, for 
example, Figure 9 of Oppenheimer & Dave 2008). The 
distance at which a wind stalls, and material begins to 
fall back, defines the size of the metal bubbles surround- 
ing individual galaxies. By circulating metals, galactic 
winds smooth out the heavy elements over the bubble. 
Multiple sightlines passing through the same bubble will 
therefore present metal-line systems that are correlated. 
The amplitude of this correlation will be roughly con- 
stant over all separations smaller than the typical size of 
the bubble. If bubbles blown by different galaxies over- 
lap, then pairs of sightlines as wide as the scale of the 
ensemble of coalesced bubbles will present a similar de- 
gree of correlation. The typical size of an enriched region 
may therefore reflect either individual bubbles or group- 
ings of bubbles. 

The size of the typical enriched region can be directly 
inferred from measurements of the spatial clustering of 
metal-line systems. On scales larger than the enriched 
regions, metal-line systems are expected to cluster like 
their sources. The correlation amplitude for galaxies de- 
pends on their halo mass and redshift, i.e., galaxy popu- 
lations with the same bias cluster the same way. Equiv- 
alently, the large-scale clustering of C IV systems deter- 
mines the bias of the galaxy population dominating en- 
richment (Scannapieco 2005). The clustering strength of 
galaxies rises towards smaller scales and steepens on sub- 
halo scales At separations smaller than the enriched 
regions, however, the clustering amplitude of metal-line 
systems will be roughly constant. The scale at which the 
C IV correlation function flattens is therefore directly re- 
lated to the typical size of the enriched regions. 

Previously, the clustering of C IV systems could only 
be measured in redshift-space along the line-of-sight. 
At large velocity separations, Scannapieco et al. (2006a) 
showed that the metal-line systems cluster like Lyman- 



9 This steepening represents the transition from the two-halo 
term to the one-halo term. It occurs on a somewhat smaller spa- 
tial scale than the break in the line-of-sight C IV auto-correlation 
function. 



break galaxies. At redshift 2.9, the clustering strength 
of Lyman-break galaxies decreases with increasing sep- 
aration as £(r) ~ (r/4.0/i _1 comoving Mpc) -16 (Adel- 
berger et al. 2005a). In contrast to this steep rise, at sep- 
arations less than 150 km s _1 , the C IV auto-correlation 
function presents nearly constant amplitude (Rauch et al. 
1996; Pichon etal. 2003; Scannapieco etal. 2006). This 
flattening characterizes measured correlation functions 
for Mg II and Si IV as well (Petitjean & Bergeron 1990; 
Scannapieco et al. 2006a) . To smooth out the metal-line 
absorption over 150 km s _1 in redshift space, enriched 
regions expanding with the Hubble flow must have phys- 
ical radii of about 500 kpc[3 The more detailed analysis 
of Scannapieco et al. (2006a) demonstrated that placing 
bubbles of metals ~ 480 physical kpc in radius around 
halos of mass ~ 10 12 M provided a reasonable fit to the 
line-of-sight correlation function. This mass scale is sim- 
ilar to the halo mass, ss 10 112 — 10 118 M Q , estimated by 
Adelberger et al. (2005a) from the clustering of Lyman- 
break galaxies. 

Such a large size for the enriched regions is surprising. 
First, it is much larger than the galactic impact param- 
eter at which strong C IV absorption is detected. The 
C IV absorption strength declines significantly beyond 
b sw 70 - 90 kpc (Steidel etal. 2010). Second, it is much 
larger than the turnaround radii of galactic winds in cos- 
mological simulations at z ~ 3 (Aguirre etal. 2001; The- 
uns etal. 2002; Kawata & Rauch 2007). In the cosmo- 
logical simulations of Oppenheimer & Dave (2008), the 
bubble radii range from 60-100 kpc over a wide redshift 
range; and the largest bubbles have not formed yet at 
z ~ 3. The analytical model of Furlanetto & Loeb (2003) 
predicts that the wind bubbles reach 100 physical kpc at 
z ~ 3. This discrepancy with the larger size indicted 
by the line-of-sight clustering of C IV systems motivates 
our measurement of C IV clustering between sightlines. 
In real space, unlike redshift space, distance and velocity 
information are not mixed. The line-of-sight correlation 
function may not directly measure physical size if en- 
riched regions have significant velocities with respect to 
the cosmic expansion. Plausible sources of these veloc- 
ity kicks include the peculiar velocities of galaxies (Davis 
etal. Peebles 1983), galactic outflows (Martin etal. 2005; 
Rupke etal. 2005; Weiner etal. 2009; Pettini etal. 2002; 
Shapley etal. 2003; Steidel etal. 2010), and cold stream 
accretion (Dekel & Birnboim 2008). 

In this paper, we provide new insight into the distribu- 
tion of carbon by comparing the redshifts of C IV systems 
toward binary quasars. Until recently, no suitable sam- 
ple of binary quasars had been compiled. Lensed sight- 
lines showed highly correlated C IV systems on interstel- 
lar scales but do not probe scales larger than a few kpc 
(e.g., Rauch etal. 1999, 2001, 2002; Ellison et al. 2004). 
At the other extreme, little correlation was seen in spec- 
tra of wide binaries with physical separations > 1.0 Mpc 
(Coppolani etal. 2006). Spectra of one binary with an 
intermediate separation of sa 50 kpc presented a number 
of C IV systems common to both spectra (Crotts et al. 
1994). Drawing targets from Hennawi (2004) and Hen- 
nawi etal. (2006a, 2006b, 2009), we obtained spectra of 
binary quasars separated by 5 to 130" thereby sampling 

10 At z ~ 3, the rate of cosmic expansion was about 
300 km s- 1 Mpc" 1 . 
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projected separations from 36 to 907 kpc, i.e., a range 
that spans the various suggested scales for the circulation 
of metals by galactic winds. Our observations differ from 
those analyzed recently by Tytler et al. (2009) in several 
important ways. Our sample has a smaller median pair 
separation (260 kpc vs about 1 Mpc). Our spectra are 
more sensitive to weak C IV systems due to their higher 
resolution. And, we examine the redshift evolution of 
the clustering properties between 1.7 < z < 4.5. (Our 
median absorption redshift is 3.0 instead of 2.0.) The 
age of the universe at redshift 4.3 is just 1.4 Myr. Sig- 
nificant enrichment and evolution may occur during the 
2.4 Gyr between redshift 4.3 and 1.7. Most importantly 
perhaps, we introduce the transverse correlation function 
to describe the clustering of C IV systems. This quanti- 
tative description of the metal distribution provides an 
empirical target (or test) for cosmological models. 

The paper is organized as follows. We describe the 
observations, data quality, and sample of C IV systems 
in § [2] In § [3l we outline several ways to count pairs 
of C IV systems, compare these pair counts to those ex- 
pected for uncorrelated absorption systems, and explain 
how the selection function of the survey was modeled. 
Appendix A derives the expressions for the error in the 
correlation function measurements. We make the criti- 
cal comparison between clustering in velocity space and 
physical space in Section [4] We discuss the size of metal- 
enriched regions and the problems with a previously used 
maximum likelihood estimator. We place an upper limit 
on the average outflow speed in § 14.21 In § 14. 3[ we con- 
strain the in-situ dispersal of metals at z ~ 3 and argue 
that the size of the enriched regions requires significant 
enrichment by winds at an earlier epoch. We propose an 
evolutionary test in § 14.41 that may empirically pin down 
the enrichment history. We summarize the conclusions 
in Section [5] discuss. 

Atomic data from Morton (2003) and a cosmology with 
Q m = 0.3, fl\ = 0.7, and H = 70 km s _1 Mpc -1 are 
used unless noted otherwise. At z — 3, these parameters 
mean that cosmic expansion produces a velocity shift of 
150 km s _1 across a physical length of 480 physical kpc, 
which corresponds to 1.35 h~ l comoving Mpc. 

2. OBSERVATIONS 

Quasars provide brighter targets for IGM tomography 
than do galaxies. The separation of lensed quasars is too 
small, however, to probe circum-galactic scales. Finding 
binary quasars with appropriate separations for CGM 
studies required systematic searches near SDSS and 2dF 
quasars (Hennawi 2004). The SDSS and 2dF quasar sam- 
ples alone contain few binaries at the relevant separa- 
tions (roughly 5 arcseconds up to an arcminute) because 
the optical fibers of spectrographs cannot be placed this 
close. At proper transverse separations / < lh^ 1 Mpc, 
Hennawi etal. (2006a) identified 221 new quasar pairs 
with 0.5 < zqso < 3.0. Another 27 high-redshift bi- 
naries, 2.9 < z < 4.3, were recently confirmed to have 
proper transverse separations from 10 kpc to 650 kpc 
(Hennawi etal. 2009). From these samples, we selected 
bright binary quasars that had similar redshifts. Our tar- 
gets are listed in Table [TJ They form 29 pairs, including 
one triplet. 

At the W. M. Keck Observatory, we obtained spectra 
of 55 quasars with the ESI spectrograph (Sheinis et al. 



2002). The echellete format provided continuous spec- 
tral coverage from 4000 A to 10000 A at 60 km s" 1 (filled 
slit) resolution. Exposure times were chosen to reach a 
S/N ratio of 15 to 20 per pixel. This spectral quality pro- 
vides sensitivity to C IV components with column density 
N(CIV) w 10 13 cm- 2 . 

The ESI spectra were reduced using the ESIRedu>0 
data reduction pipeline. Scattered light was fitted 
and subtracted from the spectroscopic flatfields. One- 
dimensional spectra were extracted for each order and av- 
eraged using a S/N weight function that varied smoothly 
across the order. Division of each spectral order by the 
fitted quasar continuum left flat, normalized orders that 
were combined into a one-dimensional spectrum. Fig- 
ure Q] shows example spectra of two quasars separated by 
100 physical kpc. 

2.1. Sample of C IV Systems 

In each quasar spectrum, we identified C IV absorp- 
tion lines by the doublet spacing. We required the two 
transitions have consistent velocity structure. In a few 
cases, this revealed blends between C IV and unrelated 
lines. We searched the bandpass between the Lya for- 
est and the observed wavelength of C IV at the quasar 
redshift. Each pair of spectra probes C IV absorption 
over a similar redshift interval because the quasars were 
selected to have roughly the same redshift. 

For illustration, the 1548 and 1551 transitions are 
marked by solid and dotted lines, respectively, in the 
J0117+3153 spectra in FigureQ] For absorption troughs 
with well separated 1548 and 155lA lines, as seen for 
example at z = 1.9987,2.0540,2.4018, and 2.4961 to- 
wards J0117+3153-A, we directly integrated the equiva- 
lent width in each line, W b s = / 1 — I(X)/I c (X)dX. We 
fitted a pair of Gaussian line profiles to measure the ab- 
sorber redshift. 

When the 1548 and 155lA lines were blended, we could 
not directly integrate their equivalent widths. Multiple 
velocity components were typically required to model 
these broad absorption troughs. We fitted Gaussian 
profiles using the SPECFIT programFi We defined 
C IV systems by linking components separated by less 
than 50 km s _1 , the effective spectral resolution (for the 
atmospheric seeing). We set the system redshift equal 
to the weighted average of the component redshifts. The 
system's equivalent width was calculated from the for- 
mula above using the fitted model of the absorption 
trough in place of the data. Linking yielded 316 inter- 
vening C IV systems and 134 associated systems. Associ- 
ated systems lie within 5000 km s _1 of the QSO (in the 
same sightline) and will not be included in the sample 
discussed in this paper. 

Figure [2] shows the redshift of intervening C IV sys- 
tems along each sightline. The paired sightlines are or- 
dered by the separation of the binary quasars. The me- 
dian sightline separation is 260 kpc with a range from 36 
kpc to 907 kpc. The absorber redshifts range from 1.7 
to 4.3 with a median of 3.0. The spectra of J1420+28, 
J1021+11, J1541+27, and J1622-I-07 provide critical cov- 
erage of small spatial scales beyond the median redshift. 

11 http://www.ucolick.org/~xavier/ESIRedux/index.html 

12 A brief description is given by G. Kriss (1994). 
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Fig. 1. — Example of sightlines separated by 100 physical kpc. 
Thick, vertical lines (red) mark Lyct and C IV at the quasar red- 
shift. Thin vertical lines (blue) mark each C IV system. We 
consider all pairs of C IV systems between sightline A and sight- 
line B. Provided Av < c, their velocity separation is Av = 
c(l-{-ZAg)~ 1 \zA — zg\, where zab is the average redshift of the two 
systems. A small velocity separation identifies coincident systems. 
The number of coincident system greatly exceeds the number we 
would find in the absence of any correlation between the sightlines. 

2.2. Distribution Function 

As discussed in detail previously (Ellison etal. 2000; 
Songaila 2001; and Scannapieco etal. 2006a), the num- 
ber of C IV systems in a quasar spectrum increases 
rapidly with improved sensitivity to column density. 
Most of the systems in our sample with rest equivalent 
width, W r (1551), less than 200 mA are optically thin. 
Their measured doublet ratio, W / & s (1548)/W / f, s (1551), 
is consistent with 2. Their C IV column density fol- 
lows directly from their equivalent width, N(C IV ) ~ 
10 140 cm- 2 (W r . i i5 5 i/200 mA). This relation provides 
only a lower limit on the column density of optically thick 
systems, however. At the modest dispersion of our spec- 
tra, optically thick lines need not be black at line cen- 
ter. We identified many optically-thick systems by their 
smaller doublet ratio, 1 < VF o6s (1548)/W / o6s (1551) < 2. 
We therefore describe the distribution of line strengths by 
their equivalent-width distribution instead of their col- 



umn density distribution. 

Distribution functions for intervening absorbers are 
usually normalized by redshift path rather than red- 
shift. A population with constant comoving space den- 
sity and constant proper size produces an equal num- 
ber of systems per unit redshift path, X(z), along a 
sightline. For the cosmology adopted here, we have 
dX = (l + z) 2 [Q,A + Cl rn (l + z) 3 }~°- 5 dz (Scannapieco et al. 
2006a) . We used our sample of intervening C IV systems 
to compute the equivalent width distribution function, 
f(W), defined as the number of systems found per unit 
redshift path per unit rest equivalent width as follows. 

For each C IV system, i, in our sample, we calcu- 
lated the total redshift path over which the system could 

have been detected, AXi — ^2j®f° AXij. Detection at 
redshift z requires the observed equivalent width of the 
weaker C IV transition exceed the 5<r measurement un- 
certainty for a system of velocity width 230 km s _1 at 
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Fig. 2. — Sightlines as a function of binary quasar separation and 
redshift probed. Each sightline is shown at the separation of the 
binary quasar. The quasars are marked by asterisks. A vertical 
tic marks C IV systems along each sightline. (Small vertical offsets 
have been applied to prevent overlapping symbols.) 

the observed wavelength^ This width is typical of a 
system near the median equivalent width in our sample. 
For each equivalent width bin, we calculate the number 
of systems per unit redshift path, 



dN 
dX 



E 



AX,, 



(1) 



where N^ is the number of systems in bin k. Normal- 
ization by the bin width, AfclogW r ,i55i = \ogW maXtk - 
logWmin,*;, yields the distribution function, 

at the mean (rest-frame) equivalent width, 

logfWOfcsNj^VlogWi, (3) 



v k £_^i - 

i=l 



13 Spectra containing the wavelength-dependence of the con- 
tinuum S/N ratio for each spectrum are available from the lead 
author. 



of systems in that bin. Figure [3J shows the resulting dis- 
tribution function. Over nearly two decades in equivalent 
width, it is well described by a power law. We attribute 
the turnover at large equivalent width, W^-,1551~ 500 niA, 
to the non-linear relation between equivalent width and 
column density. The relatively small redshift path ob- 
served at the highest sensitivity causes the large error 
bars at W r ,i55i < 20 mA. The fit to the seven points 
between 20 and 500mAis shown in Figure El This fit 
to the distribution function is used in our Monte-Carlo 
simulations in §[3] At low equivalent width, a direct com- 
parison can be made to the column density distribution 
functions. We find consistency with Ellison etal. (2000), 
Songaila (2001), and Scannapieco etal. (2006a). 

3. CORRELATION FUNCTIONS FOR C IV SYSTEMS 

Two approaches are used here to describe the cluster- 
ing of C IV systems. The first requires multiple sightlines 
at relatively small separations. The second, which has 
been applied previously, estimates separation in redshift 
space instead of real space. The two methods therefore 
measure different physical quantities. This section de- 



Martin et al. 



1000 - 



<! 100 - 







10 - 



0.1 - 



1 



100 
W r [mA] 



1000 



Fig. 3. — Equivalent width distribution function for intervening 
C IV 1551. Components have been grouped into systems using a 
linking length of 50 km s™ 1 . We fit the distribution function with 
a power law of the form, log/(A ) = b + mlog W(mA), including 
errors in W equal to one-half the bin width; see text. The uncer- 
tainty <5/k follows from the upper and lower 84.1% confidence level 
Poisson limits, corresponding to la limit for Gaussian statistics, 
on the number of objects in bin k. We find b = 4.3 ± 0.4 and 
m = -1.6 ± 0.2; or, for W in A, b = -0.6 ± 0.4. 

scribes the first method in some detail and applies both 
approaches to the data. Readers interested primarily in 
the insight gleaned from comparing the two methods may 
skip directly to § 0] 

3.1. Cross- Correlation of C IV Systems between 
Sightlines 

3.1.1. Pair Counts 

The number of absorption-line systems found at sim- 
ilar redshift, but in different sightlines, is related to the 
correlation function on the scale of the sightlines separa- 
tion. To find this relationship, data-data pairs can be de- 
fined by the apparent velocity difference of the C IV sys- 
tems in two sightlines. For example, in Figure [TJ pairs 
of C IV systems are detected in sightlines A and B at 
A 4642, A 4728, A 4950, and A 5015. Whether or not 
the system at A 5412 in sightline A and the system at A 
5430 in sightline B count as a pair depends on the ve- 
locity threshold used in the definition. In this example, 
their velocity separation is about 950 km s _1 , where 
Av = c(l + zab)~ 1 \za — zb\ and zab is the average 
redshift of the two systems. Histograms of pair counts 
are illustrated in Figure 0] for three values of the thresh- 
old velocity, Av pmr = 200,600, and 1000 km s" 1 . The 
counts of systems separated by Avab < Av pa i r have 
been summed over all pairs I of sightlines at similar trans- 
verse separation rp . For the most inclusive definition, all 
five pairs detected towards J0117+3153 contribute to the 
second bin at 300 h" 1 comoving Mpc. At fixed transverse 
separation, the number of data-data pairs declines with 
more restrictive pair definitions, i.e., smaller Av pa i r . For 
any value of Av pa i r , however, the counts in Figure 2] in- 
crease as the separation of the pairs decreases. 

We can statistically differentiate physically related 
pairs of C IV systems from chance coincidences. The 
number of unrelated pairs, hereafter random-random 
pairs RaR.b(?\p, Av pa i r ), was modeled using a Monte- 



Carlo simulation. The velocity widths of the fake sys- 
tems were fixed at 230 km s _1 . Equivalent widths were 
drawn from the distribution function shown in Figure [3) 
The redshifts of the fake systems were randomly chosen 
within the survey range for each quasar sightline. In 
the Monte Carlo simulations, fake systems were rejected 
if the observed-frame equivalent width was less than the 
5cr detection limit. Because the variance is a strong func- 
tion of observed wavelength in each spectrum, and spec- 
tral quality varied by sightline, the redshift distribution 
of the fake systems reflects the sample biases introduced 
by telluric features, quasar emission lines, and interven- 
ing absorption lines. The equivalent widths of the fake 
systems are shown in Figure [5] Their distribution ap- 
pears consistent with the real sample. Averaged over 
1000 iterations of the simulation, the average number of 
fake systems in a given sightline agrees (to within the 
Poisson error) with the actual number of C IV systems 
found. When we compare sightlines, these fake systems 
yield far fewer pairs of absorbers. 

In each bin of the histogram shown in Figure [4j the 
ratio of data-data pairs, DD(rp, Av pa i r ) to random- 
random pairs RaRb (Av pa i r ) defines a cross-correlation 
amplitude, where £,ab + 1 = DaDb/RaRb- If DaDb = 
RaRb, then all the pairs are random coincidences; and 
no clustering is detected. The pair counts are implicitly 
a function of the binning, so we drop the explicit refer- 
ence to the binning. The variance in the ratio DD/RR is 
generally larger than the Poisson error in the bin counts. 
The correction to the Poisson term is significant when 
DD > N a b s /4 (Mo etal. 1992). In most bins, the number 
of DD pairs is high compared to the number of absorbers, 
N a bs, in sightlines contributing to the bin. Calculating 
the uncertainty of the measured correlation therefore re- 
quired special attention. 

Landy & Szalay (1993) introduced the estimator, £ = 
(DD — 2DR 4- RR)/RR, to minimize the variance in a 
correlation function. For the cross-correlation of sight- 
lines, the mean number of data-random pairs is (2DR) = 
(DaRb) + (DbRa), where the average is taken over 
1000 fake sightlines with the same selection function 
as the corresponding quasar spectrum. The correlation 
strength, £ab(»"p, Av pa i r ), calculated from these pair 
counts is shown in Figure |H1 Since £ > 1 for these 
data, the Landy & Szalay estimator does not reduce the 
variance to the Poisson level. Appendix A explains how 
we computed the correction to the Poisson term. These 
larger error estimates are shown for all correlation func- 
tion figures. 

In Figure [HI the widest sightlines at rp = 2.844 h" 1 
comoving Mpc present highly correlated C IV systems. 
The systems cluster even more strongly as the sightline 
separation decreases towards 0.075 hr 1 comoving Mpc. 
The amplitude of the correlation is independent of the 
chosen binning in rp. 

As Avp a i r increases, however, the correlation becomes 
weaker in Figure [6) Stated another way, the probability 
that C IV systems (in different sightlines) are physically 
related declines as their redshift separation increases. 
This sensitivity of clustering strength to Av pa i r presents 
an obstacle for comparing the clustering of C IV systems 
to other populations. 

3.1.2. Transverse Correlation Function 
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Fig. 4. — Number of pairs of C IV systems as a function of 
sightline separation, rp. For each C IV system in sightline A, we 
computed the velocity separation between it and every C IV sys- 
tem in sightline B. Velocity separation smaller than a threshold, 



Av„ 



identify coincident systems - pairs of C IV systems at 



nearly the same rcdshift in both sightlines. We sum the num- 
ber of coincident systems over all the binary quasars I in the bin so 
that DD(rp, Av pa i r ) = J^ ; DD ! (rp, Av pa i r ). The data (open his- 
togram) present many more systems at similar redshift than would 
be expected from uncorrelated sightlines (filled histogram). The 
bin boundaries were chosen to have a similar number of random- 
random pairs, where the expectation value, RR(rp, Av pa i r ) = 

^2l(RR (rp, Avp a i r )). was obtained from a Monte Carlo simula- 
tion. From bottom to top, the three panels show how the number 
of coincident systems declines as the definition of a pair is relaxed 
to include larger velocity differences. 
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Fig. 5. — Distribution of accepted C IV W r ,i55i equivalent 
widths. Dashed, red histogram compares the fake absorbers to the 
data, represented by the solid histogram. The distribution func- 
tion for fake absorbers was cut-off at 500 mA in the Monte Carlo 
simulation, and the comparison is void for stronger systems. The 
fake and real data are consistent within their uncertainties, where 
the statistical error in the fake data is 3.1% after 1000 Monte Carlo 
iterations. 
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Fig. 6. — Cross-correlation amplitude of C IV systems vs. the 
separation of the sightlines. The value of the Landy &: Szalay 
(1993) estimator is shown at the median separation in each bin. 
The lcr error bars include a correction to the Poisson term as de- 
scribed in Appendix A. Decreasing the maximum velocity differ- 
ence defining pairs of C IV systems, A^p,^, increases the correla- 
tion amplitude between the sightlines. For each value of Av pa i r , 
the correlation amplitude declines as the sightline separation in- 
creases. 

One measure of clustering that can be evaluated for 
our sample, and compared to other populations, is the 
projection of the redshift-space correlation function. The 
projection of £,'(r) along the line-of-sight is of particular 
interest. In redshift space, the separation of two systems 
is r 2 = rp + All 2 . The distance in redshift-space is 
calculated from the velocity difference of two systems: 



AII= (l + z)Av/H(z) 



(4) 



in comoving units. The integration limits can be cho- 
sen to integrate over peculiar velocities, thereby robustly 
identifying physically related pairs of C IV systems. This 
projection 



u>(rp,II 



pair ) 



+n„ 



e(r P ,U)dU, 



(5) 



defines the transverse correlation function. It differs from 
the angular correlation function in that we use the co- 
moving distance transverse to the sightlines rp(z), rather 
than the angle between the sightlines, to describe separa- 
tion in the plane of the sky. This approach makes sense 
because Tp(z) is nearly constant over the projection in- 
terval of interest. 

Figure [7] shows the transverse correlation function. It 
is the product of ^AB{fp,/Av pa i r ) and the (comoving) 
distance 2An pa i r . Since the transverse correlation func- 
tion w has dimensions of length, we normalizing it by 
the sightlne separation and plot the dimensionless quan- 
tity w/rp. The C IV correlation amplitude increases as 
the sightline separation decreases. For the two points 
at the smallest separations, the amplitudes of w/rp arc 
essentially identical for the three projection intervals il- 
lustrated. At larger separations the correlation strength 
presents a measurable, although still small, dependence 
on the the projection interval. 

To illustrate the influence of the projection interval 
on w/rp , consider the projection of Lyman-break galax- 
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ies. Their clustering amplitude in redshift-space £'(r) is 
higher than £(r) in real-space because cosmological infall 
causes structures to appear smaller in redshift space than 
their true physical size. The reduced volume in redshift- 
space increases the clustering amplitude. To estimate 
the boost to the galaxy - galaxy correlation function, 
we apply Equation 15 of Hawkins et al. (2003) which de- 
scribes linear infall(3 The magnitude of the correction 
depends on the bias. Figure [7] shows the projection for 
b = o~8,g/o~&.cDM — 2.4, tq — 4.0 hr x comoving Mpc, and 
7 = 1.6 (Adclberger etal. 2005a). At sightline separa- 
tions > 1 ft. -1 , this infall correction noticeably boosts the 
correlation amplitude for each projection interval All; 
but it has no detectable effect at smaller separations. 
The galaxy-galaxy correlation function is not fitted to 
the C IV correlation function in Figure [7J It is therefore 
remarkable that the amount of C IV clustering agrees 
roughly with the amplitude of the galaxy - galaxy cor- 
relation function. 

Previous work revealed a strong cross-correlation am- 
plitude between Lyman-break galaxies and C IV sys- 
tems on scales > 1 h^ 1 comoving Mpc (Adelberger etal. 
2005b). We would like to explore the relationship of 
galaxies and C IV systems on the smaller scales probed 
by the binary quasars, but this turns out to be difficult 
with current data. Adelberger etal. (2005a) only mea- 
sured the galaxy - galaxy correlation function at sepa- 
rations greater than 1 hr x comoving Mpc. On sub-halo 
scales, the galaxy - galaxy correlation function is ex- 
pected to be steeper than it is in the two-halo regime. 
By extrapolating the 7 = 1.6 power-law fitted on large 
scales, the amplitude of the transverse correlation func- 
tion for Lyman-break galaxies is a lower limit at small 
separation in Figure [JJ This lower limit for galaxies 
is higher than the upper limit (one standard deviation) 
for C IV systems at small scales. In light of the uncer- 
tain magnitude of this discrepancy, this result alone is 
hardly compelling evidence for a difference in the cluster- 
ing properties of Lyman-break galaxies and C IV systems 
on small scales. 

We observed only a few binary quasars at separations 
greater than lh^ 1 comoving Mpc, so our data alone 
poorly constrain the power-law index of the correlation 
function on large scales. Taking 7 = 1.6 (like LBGs) as 
a prior, the fitted correlation length r$ = 3.65 ± 0.34 h~ l 
comoving Mpc for our sample technically agrees with the 
value r = 4.0 ± 0.6 fitted by Adelberger etal. (2005a). 
The lower value, however, reflects a best fit with a lower 
correlation amplitude overall. This difference is worth 
noting because our sightlines at separations greater than 
1 h^ 1 comoving Mpc show a slightly stronger correlation 
in Figure [7J than do galaxies. There seem to be two op- 
tions. Our measurement can be fit with a larger correla- 
tion length ro = 7.0±2.0 h^ 1 comoving Mpc and shallow 
power-law index 7 = 1.0 ± 0.2. Or, we can adopt that 
prior that C IV systems cluster like Lyman-break galax- 
ies on large scales and find that C IV do not cluster as 
strongly as Lyman-break galaxies on scales < lh^ 1 co- 
moving Mpc. To resolve this ambiguity, we re-examine 
the clustering of C IV systems along the line-of-sight be- 

14 Although linear infall theory is not strictly valid on these 
scales, Cooray & Sheth (2002) argue that it is a good approxima- 
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Fig. 7. — Transverse C IV correlation function normalized by 
sightline separation versus median sightline separation. The bins 
are identical to those in Figures f4] and [6] The lines show the pro- 
jection of the galaxy-galaxy correlation function at z ~ 2.9 (Adel- 
berger etal. 2005b), which has been extrapolated below l/i — 1 co- 
moving Mpc with constant power-law index 7 = 1.6 and perturbed 
by linear infall (Hawkins etal. 2003). Because infall reduces the 
volume in redshift space, it increases the correlation amplitude 
slightly. Independent of projection interval, however, the clustering 
amplitude of C IV systems is very similar to that of Lyman-break 
galaxies. No parameters are fitted in this figure, so their similarity 
demonstrates that these galaxies and C IV systems have nearly the 
same correlation length. 

fore returning to the best fit to the transverse correlation 
function in § 14.11 

3.2. Verification of the Line-of-Sight Correlation 
Function 

For each quasar spectrum (and fake spectrum) , we cal- 
culated the velocity separations Av of the C IV systems. 
The separations of data-data, data-random, and random- 
random pairs were binned. When we sum over all sight- 
lines I, the counts yield the correlation strength 

a ° s{Av) = E^aIu - (6) 

As discussed previously for £,ab, this estimator has a 
larger error than the Poisson term when £ > 1. Ap- 
pendix A describes our error estimate. These errors are 
shown on our correlation function in Figure [U The clus- 
tering amplitude grows steadily as the velocity difference 
of the pairs decreases from 500 to 100 km s _1 . Below 
separations of about 100 km s _1 , the clustering strength 
is constant. Our measurement confirms the break de- 
scribed by Scannapieco et al. (2006a) . Their spectra, ob- 
tained at higher spectral resolution with UVES, demon- 
strated the near constant amplitude down to the linking 
length. 

We fit the line-of-sight correlation with a simple model. 
Any model of the auto-correlation function for C IV sys- 
tems must break at small separations to fit the rela- 
tively flat correlation amplitude at small velocity sep- 
aration. Scannapieco et al. (2006a) estimated a break 
at 150 km s _1 , corresponding to 1.34 h^ 1 cMpc. We 
convert velocity differences into a redshift-path length, 
r = II = (1 + z)Av/H(z), using the median redshift of 
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Fig. 8. — Two-point correlation function of C IV absorption sys- 
tems vs. line-of-sight distance. The new ESI data (squares) agree 
with the earlier VLT data (red circles). The C IV components 
shown in Figure 10 of Scannapieco et al. (2006) were grouped into 
systems using a linking length of 50 km s" 1 . The bins near 
500 km s — 1 have a small number of systems; and we attribute the 
low pair counts to systematic error introduced by the C IV doublet 
spacing. Fixing ro = 4.0 /i —1 comoving Mpc and 7 = 1.6 in Eqn.[7] 
we fit n = 1.08 ± 0.17 h~ l comoving Mpc to the joint data (solid 
line); the dashed and dotted lines illustrate the 1-sigma error in r\. 

RR pairs in each bin. We assume a real-space correlation 
function of the form 



£(r) = (ro/r) 7 for r > r\ 
( r o/ r i) 7 for r < n. 



(7) 



The redshift-space correlation function follows from Eqn. 
15 of Hawkins etal. (2003), which estimates the pertur- 
bation caused by linear infall. Due to the slope of the cor- 
relation function, infall lowers the amplitude of the line- 
of-sight correlation function at separations much larger 
than r\ . At smaller separations where the real-space cor- 
relation amplitude is constant, the volume effect caused 
by infall again boosts the correlation amplitude (as it did 
for the transverse correlation function). 

The best fit to the joint ESI plus UVES data yields 
r = 4.9 ± 0.7 and n = 1.2 ± 0.3 for 7 = 1.6. This corre- 
lation length and power-law index are consistent with the 
values ro = 4.0±0.6/i _1 comoving Mpc and 7 = 1.6±0.1 
fitted to the galaxy-galaxy correlation function at z = 2.9 
(Adelberger etal. 2005a). Fixing r = 4.0 ft." 1 comoving 
Mpc and 7 = 1.6, we fit n = 1.08 ± 0.17 h' 1 comoving 
Mpc to the joint ESI plus UVES line-of-sight C IV cor- 
relation function. This fit and the 1-sigma uncertainty 
range are plotted in Figure [5] They provide a reasonable 
description of the data. However, the correlation length 
ro strongly influences the fitted break at r\] and we ex- 
amine the implications of this parameter degeneracy in 
the next section. 

4. DISCUSSION 

On scales larger than the typical virialized halo, the 
correlation function of galaxies is expected to represent 
the distribution of dark matter halos that host them. A 
larger correlation length ro corresponds to a greater prob- 
ability of finding a neighbor. Populations that cluster 
more strongly have higher ro and reside in more massive 



halos. From the cross-correlation of Lyman-break galax- 
ies with C IV systems, Adelberger et al. (2005b) showed 
that C IV systems with column density > 10 13 cm" 2 
cluster as strongly as Lyman-break galaxies. It follows 
that these metals and Lyman-break galaxies populate the 
same halos. Tens of thousands of Lyman-break galaxies 
have been identified photometrically, and many thousand 
observed spectroscopically. It therefore seems reasonable 
to assume that the single power law that describes the 
clustering of Lyman-break galaxies on scales greater than 
1 ft" 1 comoving Mpc also describes C IV systems on large 
scales. On scales smaller than virialized halos, multiple 
galaxies orbit within a deeper potential well; and the 
clustering of galaxies is harder to predict due to dissipa- 
tional processes. 

Our interest here is the scale of metal enriched regions 
r\. On scales r < ri, the circulation of metals by feed- 
back processes smoothes out their distribution; and the 
C IV auto-correlation amplitude is roughly constant. A 
difference between the fitted value of n in redshift space 
and in real space is of particular interest. It may reflect 
the Doppler shifts of enriched regions. 

For the C IV systems, the fitted v\ depends on the cor- 
relation length. As r increases, the clustering amplitude 
at any separation r > ro increases. To remain consistent 
with the measured correlation strength at small separa- 
tions, the n transition must shift to larger scales, where 
the correlation amplitude is lower. Figure |H] shows how 
the fitted t\ grows with increasing correlation length. 
Based on the shape of the contours, the line-of-sight mea- 
surements clearly constrain the scale r± more tightly than 
they do the correlation length ro. The difference in r\ 
values and their implied Doppler shifts decrease as the 
correlation length increases. 

The underlying real-space correlation length ro must 
describe both the transverse and line-of-sight clustering 
measurements. Clustering along the line-of-sight pro- 
vides the better constraint because our sightlines to the 
binary quasars were simply not far enough apart to con- 
strain ro well. Both measured C IV correlation func- 
tions, nonetheless, are marginally consistent with the 
Lyman-break galaxy correlation length, vertical line at 
ro = 4.0 ± 0.6 ft" 1 comoving Mpc in Figure |U The best- 
fit correlation length ro = 4.9±0.7ft _1 is slightly larger. 
In the absence of any peculiar velocities, the correlation 
length must be increased to 6 ft -1 comoving Mpc before 
the same t\ fits the clustering measured in both redshift- 
space and real space. Our C IV observations provide 
only weak constraints on very large-scale clustering, so 
they alone do not rule out ro ~ 6 ft -1 comoving Mpc. 
However, such halos are about three times more massive 
than those surrounding Lyman-break galaxies and have 
a number density about 4 times lower Fi In light of the 
strong cross-correlation amplitude between C IV systems 
and galaxies (Adelberger et al. 2005b) and the presence of 
outflows from most Lyman-break galaxies (Steidel et al. 
2010), we do not discuss this possibility further. 

For ro = 4.0 or 4.9 ft -1 comoving Mpc, the implied pe- 
culiar velocities set an upper limit on the average speed 
of galactic outflows. Together, the size and speed mea- 
surements place interesting new constraints on the era of 
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Fig. 9. — Contours of Ax 2 for two-parameter fits to the C IV cor- 
relation function. The correlation length ro and the size of enriched 
regions r\ parameterize Eqn.0 and the power-law index is taken to 
be 7 = 1.6. The fitted transverse and line-of-sight C IV correlation 
functions are illustrated separately in the upper and lower panels, 
respectively. The contours indicate values of Ax 2 , where 2.30, 4.61, 
and 6.17 correspond to the 68.3%, 90.0%, and 95.4% confidence 
level. The line-of-sight and transverse correlation functions sam- 
ple the same enriched regions, so they must be jointly fit with the 
same ro and r\ . Adopting the clustering length of Lyman-break 
galaxies (vertical line) implies the presence of peculiar velocities 
among C IV systems. Larger values of the correlation length re- 
quire smaller velocities. No Doppler shifts are required when the 
correlation length increases to (the rather large value of) 6/i — 1 co- 
moving Mpc and r\ ~ 1.3 ft -1 comoving Mpc. (The fits to the 
line-of-sight clustering bifurcate at ro>4.5/i -1 comoving Mpc due 
to errors in counting C IV systems separated by the doublet spac- 
ing, Ad = 500 km s~ x .) 

metal dispersal. 

4.1. The Size of Metal-Enriched Regions 
4.1.1. Amplitude of the Transverse Correlation Function 

As our best estimate for the size of metal-enriched re- 
gions, we adopt ro = 4.0 ft. -1 comoving Mpc and 7 = 1.6 
(Adelberger etal. 2005a). We fit the model given by 
Equation[7]to the transverse correlation function and find 
the power-law flattens on a scale n = 0.42 ± 0.15 ft -1 
comoving Mpc. This fit is shown in the top panel of 
Figure QUI where it provides a good description of our 
measurements. For ro = 4.0 ft. -1 comoving Mpc, the fit 
to the line-of-sight correlation function shown in Figure |8] 
requires larger v\ = 1.08 ± 0.17 ft -1 comoving Mpc. A 
model illustrating this larger transition scale is shown in 
the bottom panel of Figure I1Q[ where it clearly under- 
predicts the correlation between sightlines seen on small 
scales. We find C IV systems at similar redshifts in the 
two sightlines more often than we would in a model with 
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Fig. 10. — Model of the correlation function defined by Equa- 
tion[7]in real-space and then corrected for linear infall using Equa- 
tion 15 of Hawkins etal. (2003). We fix the correlation length 
at ro = 4.0 ft -1 comoving Mpc and power-law index at 7 = 1.6, 
(Adelberger et al. 2005a). We illustrate how the scale n affects the 
transverse correlation amplitude. Data are as shown in Fig. [7] (a, 
Top) A fit to these data (specifically Av = 100 km s" 1 ) suggests 
ri = 0.42 ± 0.15 ft" 1 comoving Mpc. We show this model and in- 
terpret ri as the typical radius of an enriched region, (b, Bottom) 
We show a model with n = 1.08 /i —1 comoving Mpc, as favored 
by our fit to the line-of-sight correlation function. The larger size 
for the enriched regions provide too little power on small scales. 

metal bubbles as large as those suggested by Scanna- 
pieco etal. (2006a). Consistency requires a model with 
redshift-space distortions. 

We also illustrate the discrepancy in redshift-space 
when the C IV correlation length takes our best-fit value 
of ro = 4.9 ft -1 comoving Mpc for 7 = 1.6. We fit- 
ted r\ = 1.21 ± 0.25 ft. -1 comoving Mpc to the line-of- 
sight correlation function. To the transverse measure- 
ment, we fitted n = 0.74 ±0.21 ft -1 comoving Mpc. The 
scale of the line-of-sight break remains significantly larger 
than the fitted transverse break. Evidently, the metal- 
enriched regions present peculiar velocities with respect 
to the cosmic expansion. We fit the magnitude of these 
velocity kicks in § 14.21 

4.1.2. Coherence Length from Coincident Absorption in 
Spectra of Binary Quasars 

The spectra of binary or lensed quasars have been used 
previously to constrain the transverse size of absorbing 
structures. The transverse size of a coherent region was 
inferred to be at least as large as the maximum sepa- 
ration of common spectral features (McGill 1990; Lopez 
etal. 2000; Ellison etal. 2004). As implemented previ- 
ously, however, the resulting coherence length has little 
to do with the size of enriched regions. To illustrate 
this important point, we apply the standard coherence- 
length analysis to our sample of C IV systems, explain 
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why it fails, and suggest a revised technique. This re- 
vised coherence-length estimator can be compared to the 
size scale r\ — 0.42 or 0.74 h~ l comoving Mpc fitted to 
the transverse correlation function. 

To estimate the coherence scale of C IV absorption, 
we first follow the formalism of McGill (1990) with the 
modification suggested by Dinshaw etal. (1997). Their 
argument is entirely geometric in nature. They describe 
absorption-line systems with a non-evolving population 
of clouds. For illustration, we adopt a spherical distribu- 
tion of clouds, but the method easily generalizes to cylin- 
drical filaments or inclined disks. The radius R describes 
the coherence scale of the cloud distribution. For pairs 
of sightlines separated by a distance Z0 the probability 
that the second sightline intersects the sphere, given that 
the first sightline "hit" the sphere is 



Pi = 2/ttIcos- 1 X(z) - X(z)yjl - X(z) 2 }, (8) 

where X(z) = l(z)/(2R) and < X(z) < 1. If either 
sightline shows absorption, the probability of coincident 
absorption, defined as C IV absorption in both sightlines 
at similar redshift, is 



P 2 = Pi/(2 - P x ). 



(9) 



The probability of absorption in just one sightline, an 
anti-coincidence, is 1 — P 2 . For a particular coherence 
length _R, the likelihood of obtaining the observed num- 
ber of coincident systems and the observed number of 
anti-coincidences is 



r N, 



c(R) = nfj 1 p 2 [x^ i )]n;i 1 (i - p 2 [x( Zj )]), (10) 

where the two products are over all coincidences Nt 
and anti-coincidences Nj, respectively. Figure [TT1 shows 
the likelihood for many different values of the coherence 
length. 

Applied to the entire sample of 29 binary quasars, the 
maximum likelihood is found at the very large scale of 
R = 1.07JIq'o 2 /i _1 comoving Mpc. Detection of just a 
single coincidence between sightlines separated by I = 
2.1 h^ 1 comoving Mpc sets a sharp lower bound. In gen- 
eral, the analysis proposed by McGill (1990) and Din- 
shaw etal. (1997) requires the sphere diameter to be at 
least as large as the separation of the widest pair of sight- 
lines with a coincidence. As we argued at the beginning 
of §|H however, the sources of metals are highly clustered. 
When two sightlines are separated by ~ 1 h" 1 comoving 
Mpc, the correlation between the halos of their sources 
causes some absorption at a common redshift. The clus- 
tering of galaxies, and therefore metal-enriched regions, 
has not been accounted for in Eqn. [8] This oversight 
leads to a coherence scale that grows with the separation 
of the widest binary quasar in the sample. Figure [TT] 
demonstrates this lack of convergence by comparing the 
likelihood functions for three subsets of data. 

In principle, it might be possible to modify the analysis 
to obtain a robust coherence length for enriched regions. 
This provides no obvious advantage over the C IV corre- 
lation function but offers a sanity check on the inferences 

16 Although previous authors used the physical separation of 
the sightlines, the co-moving separation must be used when the 
absorber redshifts span a significant range in the growth factor, 

1 + z. 



drawn from the correlation function. We suggest that 
the maximum-likelihood estimate of the coherence length 
may converge provided the sightline separations remain 
smaller than the typical size of individual enriched re- 
gions. When the likelihood analysis yields R compara- 
ble to one-half the maximum sightline separation, the 
method has not converged. To estimate this coherence 
scale, we calculated the function given by Eqn. [10] for a 
series of 1, 2, 3, ..., and 29 binary quasars. Figure [T2l 
compares our result to the line 2R — I. When the sight- 
line separation is smaller than the diameter of the typical 
enriched region, / < 2R, we expect many coincidences, 
few anti-coincidences, and a most-likely coherence diam- 
eter 2R significantly larger than the maximum sightline 
separation. 

In Figure [12l the coherence length estimate for all sub- 
samples with maximum sightline separation less than 
0.6 ft, -1 comoving Mpc allows R « 0.3 hr 1 comoving 
Mpc. In contrast, as the maximum sightline separation 
grows larger than 0.6 ft -1 comoving Mpc, the estimated 
coherence length follows the locus R — 0.5Z. The coher- 
ence length appears to converge to a value R « 0.3 hr 1 
comoving Mpc in Figure Q21 This scale is consistent with 
the 0.42±0.15 hr 1 comoving Mpc fitted to the transverse 
correlation function with r§ — 4.0 h~ x comoving Mpc. It 
is marginally inconsistent with the larger n fitted when 
a larger correlation length is adopted. Given the many 
implicit assumptions in the interpretation of both R and 
r\ , and the acknowledgement that enriched regions could 
have a range of sizes, we believe the revised coherence 
length analysis lends strong support to our proposed in- 
terpretation of the C IV correlation function. It may rule 
out n as large as 0.7 h^ 1 comoving Mpc, but we prefer to 
see the technique calibrated against cosmological models 
before making this claim. 

4.1.3. Previous Observations of C IV Absorption near 
Lyman-Break Galaxies 

Sightlines passing by known Lyman-break galaxies di- 
rectly probe the properties of circum-galactic gas. Spec- 
tra of background light beacons present intervening ab- 
sorption near the redshift of Lyman-break galaxies. Be- 
cause background galaxies are numerous, spectroscopy 
of galaxies can probe circum-galactic gas at low impact 
parameter. Background quasars are less common, but 
spectra of the brightest quasars provide higher sensitiv- 
ity to absorption lines. The impact parameters of these 
detections can be compared to our measured scale for the 
enriched regions, about 150 physical kpc. 

Composite galaxy spectra were recently created over a 
range of impact parameters, b, from Lyman-break galax- 
ies (Steidel etal. 2010). Intervening C IV absorption was 
detected at b = 60 physical kpc but not at b = 100 kpc. 
Our analysis suggests the radius of the typical metal- 
enriched region is significantly larger. Is this a sensitivity 
difference or a discrepancy? 

Our sample contains C IV systems as strong as the 
Wr,i549 = 700 mA absorption trough in the composite 
spectrum at b — 60 kpc. Our spectra of median qual- 
ity are more sensitive than their detection threshold of 
Wr,i549 = 120 mA at b — 100 kpc. Figure |3] illustrates 
our equivalent width distribution. The typical detection 
limit is W r i55i = 0.5T4 7 r ,i549 ~ 20 mA ; and a few spec- 
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Fig. 11. — Likelihood of observed coincident and anti-coincident 
C IV systems (normalized by the maximum likelihood) vs. the co- 
herence length of a spherical distribution of clouds. The dashed 
line shows the cumlative distribution, which we use to define the 
95% confidence limit on the statistical error. We defined coinci- 
dences as systems in different sightlines within 600 km s _1 of each 



other. The most likely coherence radius R 



1.07"' 



, h,- 1 Mpc 



for all 29 pairs. However, the best radius drops to R ~ 0.5 h 
Mpc when the 9 widest binaries are dropped from the sample; and 
R ~ 0.3 h —1 Mpc when the analysis is applied to only the clos- 
est 8 binaries in the sample. As sample size increases (and wider 
binaries are included), the maximum-likelihood analysis does not 
converge to a coherence length. See text for explanation. 
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Fig. 12. — Coherence length (in comoving units) for different 
subsets of paired sightlines. Using our large set of pair observations 
at varying separations, we iteratively threw out the largest pair and 
recomputed the coherence length using the maximum likelihood 
analysis. The x — axis inidicates the widest (comoving) separation 
included in the subsample. The diagonal line indicates a sphere 
radius of one-half the maximum sightline separation, R = 0.51. The 
estimated coherence length does not converge as the sample size 
is increased. When R is significantly greater than 0.5Z, however, 
this method may provide a robust estimate of the coherence length 
(shaded, red) of individual enriched regions. 

tra reach W r ^ 5 i = 0.5W r . ! i 5 4g ~ 8 mA . Our sensitivity 
to weaker systems may explain why we find a larger size 
for the enriched regions. 

These weaker C IV systems do appear to be associated 
with the environment of Lyman-break galaxies. The clus- 



tering of weaker systems has been measured using sen- 
sitive spectra of 23 quasars in Lyman-break galaxy sur- 
vey fields. The cross-correlation length (between LBGs 
and C IV systems) increases slowly with the C IV column 
density (Adelberger et al. 2005b, Figure 11). At the com- 
pleteness limit of our survey, W^issi ~ 20 mA, C IV sys- 
tems cluster as strongly with Lyman-break galaxies as 
Lyman-break galaxies do with each other. Hence, we can 
be reasonably confident that the large metal enriched re- 
gions describe the environment of Lyman-break galaxies. 
We have not shown, nor do we claim, that these weaker 
C IV systems mark outflows from Lyman-break galaxies. 

4.2. The Peculiar Velocities of Metal-Enriched Regions 

It is clearly of interest to understand why the line- 
of-sight correlation between C IV systems flattens on a 
larger scale than does their transverse correlation func- 
tion. The line-of-sight correlation function is derived 
from the redshift difference of pairs of C IV systems. 
Their separation in redshift mixes information about 
their physical separation and their relative Doppler 
shifts. Non-zero peculiar velocities smear out the cor- 
relation in redshift space. Adding random velocity kicks 
to enriched regions therefore flattens the line-of-sight cor- 
relation function on a scale somewhat larger than r\ . We 
can fit the magnitude of the typical velocity kick. 

As suggested by Hawkins etal. (2003), we represent 
the random motions by an exponential of the form 



f(v) 



1 



aV% 



exp 



V2\v\ 



(11) 



The impact of the Doppler velocity a on the line-of-sight 
correlation function is shown in Figure 1131 where these 
kicks increase from zero to a = 1000 km s _1 . The 
correlation amplitude decreases at velocity separations 
Av < a and grows on the scale Av ~ a. The magni- 
tude of the velocity kick required to fit the correlation 
function depends on j-j., rn, and 7. The largest value 
of a results if we adopt the smallest value of r\ (and 
therefore the smallest correlation length). For example, 
taking r\ = 0.42 h^ 1 comoving Mpc as shown in the top 
panel of FigurcQ21 we fitted a ~ 120 km s _1 - the model 
shown by the bold line. With no velocity kick, wc fitted 
a larger n = 1.08 hr 1 comoving Mpc to these data in 
Figured] Models with these larger enriched regions do 
not allow any velocity kicks. As illustrated in the bottom 
panel of Figure 1131 large enriched regions with peculiar 
velocities yield too few pairs of C IV systems at sepa- 
rations of a few hundred km s _1 . To summarize, we 
quantify the upper limit on a in Figure 1141 Varying the 
Doppler velocity and computing the fit statistic, we find 
a 3<7 upper limit a < 300 km s _1 . We conclude that the 
average Doppler shift of metal- enriched regions along the 
line-of-sight is likely about 120 km s _1 and appears very 
unlikely to exceed 300 km s _1 . 

The average peculiar velocity of metal-enriched regions 
in three-dimensions is v pec — \/3a ~ 200 km s _1 (3er up- 
per limit of 500 km s _1 ) after correction for projection 
along the line-of-sight. These peculiar velocities are sim- 
ilar to those of galaxies. Their primary source may be 
gravitational instability associated with the underlying 
mass structure. The outflow speed measured in galaxy 
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spectra is ~ 200 km s -1 (Martin 2005; Weiner et al. 2009; 
Steidel etal. 2010). This absorbing gas lies close the 
host galaxy (Martin & Bouche 2009). The average out- 
flow speed, where the average is a spatial average from 
the starburst region out to the bubble radius where the 
wind stalls or even turns around, has not been directly 
measured previously. Our measurement shows that this 
average outflow speed is very unlikely to be larger than 
200 km s -1 . The average velocity may be consider- 
ably less because gravitationally-induced peculiar veloc- 
ities must contribute to random velocity kicks. 

Adelberger et al. (2005b) reached a very different con- 
clusion about average outflow speeds. Their Figure 12 
shows the distribution of redshift differences between 
C IV systems and Lyman-break galaxies. If the galaxy 
and gas are co-spatial, the typical offset in redshift corre- 
sponds to a Dopplcr shift of ~ 500 km s -1 . Because such 
velocities are larger than virial motions within Lyman- 
break galaxies, this observation was interpreted as evi- 
dence for galactic outflows at 500 km s -1 or more. 

Outflow speeds of 500 km s _1 are inconsistent with the 
measured line-of-sight correlation function for C IV sys- 
tems. We suggest that the redshift differences between 
C IV systems and Lyman-break galaxies can be under- 
stood without appealing to fast outflows. The physi- 
cal separation between the absorbing gas and the galaxy 
may explain their redshift differences. For example, cos- 
mic expansion at z ~ 3 produces a velocity difference of 
450 km s" 1 over the 4.0 ft. -1 comoving Mpc (the corre- 
lation length). As the column density of C IV systems 
increases, their cross-correlation length (with galaxies) 
increases. For larger cross-correlation length, which is 
equivalent to a stronger correlation amplitude, the dis- 
tance from a galaxy to the nearest intervening systems 
is smaller. The average column density of the absorp- 
tion will be higher the smaller the separation or impact 
parameter. The strongest C IV sy terns are always found 
at small impact parameters, so random spatial offsets 
between galaxies and C IV systems can plausibly explain 
the increase in redshift differences observed with growing 
impact parameter. 

4.3. The Dispersal of Metals 

Understanding how C IV systems cluster reshapes our 
view of the IGM enrichment history. In particular, the 
large sizes and low peculiar velocities of enriched regions 
constrain the era of metal dispersal. This process is of 
course on-going as there is strong observational evidence 
for galactic winds across much of cosmic time. Not all 
galaxy populations, however, can easily distribute metals 
as widely as observed. 

4.3.1. Lyman-break Galaxies at z ~ 3 

How far can a galactic wind from an individual galaxy 
can transport metals? In the absence of any forces to 
stall the expansion, the maximum distance traveled is 
limited by the age of the starburst galaxy and average ve- 
locity. Adopting the upper limit (v) ~ 200 km s -1 com- 
patible with line-of-sight clustering observations, galactic 
winds travel 
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Fig. 13. — Models for t he line-of-sight correlation function. The 
parameter a in Eqn. 1111 describes the typical random velocity of 
an enriched region. The real-space correlat ion function and infall 
models are the same as those shown in Fig. 1101 (with ro = 4.0 /i -1 
comoving Mpc and 7 = 1.6). Data are as in Figure \E\ with the 
bins affected by the doublet spacing omitted for clarity. (Top) 
We fix r\ = 0.42 h~ 1 comoving Mpc, the value favored by our fit 
to the transverse correlation function. Thick line shows the best 
fit peculiar velocity a = 120 km s — 1 . Thin lines illustrate other 
peculiar velocities. Lower velocities give too much correlation on 
small scales; and higher velocities require more pairs at velocity 
separations An ss 300 to 1000 km s -1 than are observed. (Bottom) 
Here we fix r\ = 1.08 h~ 1 comoving Mpc, the fitted bubble size for 
no peculiar velocity (bold line). We show that adding even a small 
peculiar velocity of a = 120 km s~ x results in too few pairs of 
C IV systems at Av < a and too many pairs at Av > a. 

in a typical LBG lifetime^ This distance is at most 
40% of our estimated radius for an enriched region. 
While typical Lyman-break galaxies clearly eject metals 
in winds, other galaxies may spread metals over larger 
regions of intergalactic space. 
The properties of Lyman-break galaxies do span a sign- 

17 For purposes of illustration, we use the median timescale for 
star formation, 321 Myr, from Table 4 of Shapley etal. (2001). 
With only optical photometry, the absolute ages are quite sensitive 
to the model for the star formation history. For example, Sawicki & 
Yee (1998) fit much younger ages ~ 30 Myr. As shown by Shapley 
etal. (2001), measurement of the Balmer break distinguishes these 
interpretations and favors typical ages of a few hundred Myr. 
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Fig. 14. — Values of peculiar velocity a fitted to the line-of-sight 
C IV correlation function. With n fixed at 0.42 h _1 comoving Mpc 
(solid line), the deviation in x 2 from the best fit allows peculiar 
velocities up to 170, 220, and 260 km s — 1 at the 1, 2, and 3<x 
confidence level. A more general 2-parameter fit to a and rj si- 
multaneously (dashed line) allows a maximum peculiar velocity of 
310 km s — 1 at the 3a level. (The correlation function of Adelberger 
etal. (2005a) is assumed for ro and 7.) 

ficiant range; and the winds from those with the largest 
lifetimes can enrich the large regions we find. To disperse 
metals out to 200 kpc, the required age is ~ 1 Gyr for 
an average speed ~ 200 km s _1 . Ages ~ 1 Gyr are fit to 
about 20% of Lyman-break galaxies (Shapley et al. 2001). 
Interestingly, the Lyman-break galaxies with older ages 
tend to be less luminous (Shapley etal. 2001). It follows 
that the lowest luminosity Lyman-break galaxies (in cur- 
rent samples) could be the source of the most dispersed 
metals. 

Given the absorption cross-section of bubbles with 
ri ~ 0.42 hr x comoving Mpc, the Lyman-break galax- 
ies are not numerous enough to be the source of most 
of the C IV systems. The number density of Lyman- 
break galaxies is n g sa 3.7 x 10 -3 h 3 Mpc -3 (Steidel 
et al. 2010). If we could place these bubbles around every 
Lyman-break galaxy, the resulting redshift-path density 
of C IV systems would be about 



dN 
dz 
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0.42 h- 1 Mpc 
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For comparison, integration of our distribution function 
shown in Figure [3] to W r ,i55i ~ 20 mA yields about 14, 
C IV systems per unit redshift. Each bubble of metals 
might yield more than one C IV system, but C IV sys- 
tems outnumber Lyman-break galaxies by an order of 
magnitude. 

4.3.2. Post-Reionization Dispersal of Metals 

Large bubbles of metals surrounding Lyman-break 
galaxies could be relics of an era when galaxies were 
smaller, more numerous, and closer together (Porciani 
& Madau 2005; Scannapieco 2005). How do the new 
results on the clustering, size, and peculiar velocities of 



enriched regions fit into this picture? 

In order for their remnant winds to cluster like 
C IV systems, higher-redshift sources must have lower 
mass. Consider the post-reionization galaxy population 
for example. Redshift 6 galaxies in halos about 30 times 
less massive than those of Lyman-break galaxies have the 
same bias as the Lyman-break galaxies at z ~ 3. At red- 
shift 4.3 galaxies 5 times less massive than Lyman-break 
galaxies would eject metals that clustered like Lyman- 
break galaxies by z ~ 3. 

Low-mass galaxies naturally produce large enriched re- 
gions. In a shallow gravitational potential, gas retain- 
ment is more difficult. Because enrichment is earlier, 
the bubbles around low-mass galaxies have more time 
to grow. At an average speed of ~ 200 km s _1 , the 
bubble could travel 150 physical kpc between redshift 
4.3 and 3. If growth stalls relative to the surrounding 
medium, the bubbles continue to grow in physical size 
due to cosmic expansion. (For a given physical size, the 
comoving bubble scale is larger the earlier the enrich- 
ment.) Most importantly, however, the number density 
of halos (and presumably galaxies) grows rapidly with de- 
creasing mass. Those low-mass galaxies between z = 4.3 
and 6 that cluster like Lyman-break galaxies could easily 
be 5 to 30 times more numerous. 

In the post-reionization universe, the potentially high 
number of galactic winds and large cross-section for 
absorption yield a large density of metal-line systems. 
From Eqn. 1131 the number per unit redshift is about 5 
when starburst galaxies at redshift 4.3 are the source 
of metals observed at redshift 3. The redshift-path 
density grows rapidly for earlier dispersal. The spatial 
distribution of the entire population of C IV systems 
with N(C IV )> 10 13 cm~ 2 can easily be attributed to 
their relic galactic winds. Significant metal enrichment 
will of course continue until much later, but the post- 
reionization galaxies appear to dictate the size scale of 
the circumgalactic medium. 

The bubbles from numerous low mass galaxies may 
overlap. Bubble coalescence would generate rapid growth 
in the size of the enriched regions. Such growth is demon- 
strated by simple analytic models. In Figure 2 of Scan- 
napieco etal. (2005), the comoving size of enriched re- 
gions grows slowly when rare galaxies provide the met- 
als. Indeed, the bubbles blown by different Lyman- 
break galaxies never overlap. The probability of find- 
ing a second galaxy within the bubble volume V is 
P = n ga iV(l + £( r i))- Substituting typical values, we 
estimate 



P = 0.04 



ngal 

3.7 x 10- 3 /i 3 Mpc" 

\V l + g(ri) 

1 + 37 



0.42 h- 1 Mpc 



(14) 



where we adopt a correlation length of 4.0 h~ l comoving 
Mpc and 7 = 1.6 to estimate the probability of finding 
another Lyman-break galaxy within 7*1 h~ l of the first. 
Since the bubbles blown by a Lyman-break galaxy with 
median age are not this large, this parameterization gives 
an upper limit on the odds of overlap. The z ~ 4.3 star- 
bursts are the youngest population (750 Myr at z ~ 3) 
that have time to blow the large bubbles we found. Their 
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5 times higher source density yields P«0.2 suggesting 
that many bubbles might remain isolated although ap- 
proaching overlap. In contrast, relic winds from redshift 

6 sources would overlap create enriched regions much 
larger than the bubbles blown by individual galaxies. 

In summary, the measured size r\ of the enriched re- 
gions is difficult to explain by metal dispersal within just 
a few hundred Myr of z ~ 3. This observation does not 
distinguish slightly earlier enrichment at z ~ 4.3 from en- 
richment at the end of cosmic reionization at z k, 6. En- 
richment over this era, however, might be distinguished 
from the Pop III scenario by measuring rapid growth in 
the size of enriched regions. Pop III stars in 1000 K ha- 
los would produce rapid growth in the comoving size of 
the enriched regions prior to redshift 6. Post- reionization 
wind bubbles, on the other hand, may coalesce between 
redshift 6 and 3, generating rapid growth in the size of en- 
riched regions. Our measured size, although large, does 
not require coalescence. Starbursts at redshifts 4-5 could 
create such bubbles, which would not have merged for the 
most part by z ~ 3. Some of the parent galaxies would 
still be observable at redshift 3 and could be identified 
by their older stellar population and low mass. 

4.4. Evolution of the Size of Enriched Regions 

The evolution of the C IV auto-correlation function 
may provide a useful constraint on the metal enrichment 
history. The size of enriched regions can only grow in co- 
moving units if bubbles merge or additional enrichment 
takes place. As bubbles coalesce due to galactic infall, 
the size of the enriched region grows as gravity assembles 
galaxies into large-scale structures. The correlation am- 
plitude for galaxy halos grows as (1 + z)~ 2 . Evolution in 
the C IV correlation function inconsistent with (1 + z)~ 2 
would indicate the addition or removal of metals. 

Scannapieco et al. (2006a) found that the amplitude 
of the line-of-sight correlation functions for Mg II and 
Fe II , and possibly C IV , increases with cosmic time. 
We used our median redshift to define lower redshift 
C IV systems (1.7 < z a b < 3) and higher redshift sys- 
tems (3 < z a b < 4.3). For each subsample, we counted 
DaDa , DaRa j and RaRa pairs in a single velocity bin 
at 50 < Av < 500 km s _1 . The median pair separa- 
tion is about Av sw 150 km s _1 at both redshifts. The 
correlation amplitude of the higher redshift subsample is 
£ + 1 ~ 5.5. The amplitude grows to £ + 1 sw 11 for lower 
redshift subsample. Hence, we confirm the amplitude of 
£los(CIV), increases with cosmic time. 

When we compare two sightlines, we count pairs with 
separations up to Av pa i r = 1000 km s _1 . The identi- 
fication of systems in different spectra is not biased by 
the C IV doublet spacing, and the larger velocity differ- 
ence reduces the noise. In the higher redshift subsample, 
the median redshift of an RaRb pair is 3.33. The me- 
dian pair redshift is 2.46 in the lower redshift subsam- 
ple. We calculate the transverse correlation function as 
described in § 13.1.21 and show the result in Figure [TBI 
At all sightline separations, the amplitude for the lower- 
redshift C IV systems is larger than that for the higher- 
redshift systems. The magnitude of the increase in clus- 
tering strength is consistent with that measured along 
the line-of-sight. We estimate £(2.46)/£(3.33) « 2 for 
C IV systems. 

Between the median subsample redshifts of 3.33 and 



2.46, we expect the amplitude of of the correlation to 
grow by a factor 1.6. Our measured correlation functions 
for the subsamples are consistent with this growth rate, 
although the large error bars allow somewhat faster or 
slower growth. Measuring evolution in the scale of the 
break r% may be possible in the future. With the current 
sample, the uncertainties in the fitted r\ are large and 
consistent with no evolution in r% larger than a factor 
of two. The addition of larger separation pairs to our 
analysis is important to allow fitting r$ and n without 
additional priors. 
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Fig. 15. — Transverse correlation function, as in Figure [7] with 
Av < 1000 km s — 1 , but for high and low redshift C IV systems. 
The correlation amplitude for the lower redshift subsample at z = 
2.46 (blue, half-filled symbols) exceeds that measured for higher 
redshift subsample at z = 3.33 (red, starred symbols). 



5. SUMMARY 

We obtained echellete spectra of binary quasars, which 
probe the IGM along parallel sightlines. Metal-enriched 
regions were identified by C IV absorption systems. We 
measured their clustering in real space, transverse to the 
sightlines, and in redshift space, along individual sight- 
lines. At large separations, the clustering of interven- 
ing C IV systems is consistent with the spatial distri- 
bution of galaxies. (Galaxies, as used here, referred to 
Lyman-break galaxies at z ~ 2.9 with correlation length 
ro = 4.0 h^ 1 comoving Mpc.) We found metals more 
smoothly distributed than galaxies, however, below a 
scale that we define as the enriched region. This scale 
may represent the turnaround radii of galactic winds or 
the overlap of multiple bubbles. We showed that fast 
outflows would smear out the C IV correlation function 
in redshift space. By comparing the correlation func- 
tions in redshift space and real space, we found the maxi- 
mum root-mean-square outflow velocity. We argued that 
the large size of the enriched regions and slow outflow 
speeds are incompatible with dispersal by winds from 
z ~ Lyman-break galaxies. 

We list our empirical results here before discussing 
their implications. 

• We measured redshifts and equivalent widths for 
316 intervening C IV systems at 1.7 < z < 4.3. 
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The number of systems per unit redshift path falls 
as the strength of the absorption increases. The 
distribution of equivalent widths is consistent with 
previous measurements of the column density dis- 
tribution (Ellison etal. 2000; Songaila 2001; and 
Scannapieco etal. 2006a). 

• We counted pairs of C IV systems as a function 
of their separation. The ratio of these counts to 
the number of pairs in simulated sightlincs deter- 
mined the clustering strength (at each separation). 
The large clustering amplitude, £ > 1, required 
special treatment of the errors. In an Appendix, 
we outlined how to compute a correction to the 
Poisson term that characterizes the Landy & Sza- 
lay (1993) correlation estimator. We learned that 
reducing the uncertainty at fixed separation will 
require more systems per sightlinc, i.e., more sen- 
sitive observations rather than more quasars. 

• We examined the correlation of C IV systems be- 
tween sightlines. When a C IV system is present in 
sightline A, we found a better than average chance 
of finding a C IV system within 1000 km s _1 in 
sightline B. The correlation amplitude grows as 
the separation of the sightlines decreases. We 
presented a projection of the C IV correlation 
function over transverse scales from 0.075 < 
r p (h~ 1 comoving Mpc) < 2.84. The fitted break 
scale, which defines the radius of an enriched re- 
gion, was found to be sensitive to the correla- 
tion length. Taking the prior tq = 4.0 hr x co- 
moving Mpc, we found r\ — 0.42 ± 0.15 ft -1 co- 
moving Mpc. However, larger enriched regions 
and correlation lengths improved the fit statistic 
slightly, e.g., r\ = 0.74 ± 0.21 hr 1 comoving Mpc 
and ro = 4.9 h~ x comoving Mpc. Repeating our 
analysis with more widely separated sightlines will 
break this degeneracy. 

• The traditional method for estimating the coher- 
ence length of enriched regions breaks down due to 
galaxy-galaxy clustering. We demonstrated that 
the result cannot be trusted when the most likely 
diameter equals the separation of the widest bi- 
nary presenting C IV systems at similar redshift. 
We propose a revised maximum-likelihood analy- 
sis. When applied to our data, it converges to a 
coherence length of 0.3 h~ l comoving Mpc. This 
result is consistent with our interpretation of the 
transverse correlation function. 

• We measured the clustering of C IV systems along 
single sightlines. We confirmed that the correla- 
tion function flattens at small velocity separation 
as suggested previously (Rauch et al. 1996; Pichon 
etal. 2003; Scannapieco etal. 2006a). Neglecting 
peculiar velocities, we fitted a correlation length, 
ro = 4.9 ± 0.7 /i _1 comoving Mpc, and break, 
n = 1.2 ± 0.3 ft- -1 comoving Mpc (where 7 = 1.6). 
Since the correlation length was consistent with 
that fitted to galaxies, we adopted tq = 4.0 hr 1 
comoving Mpc (Adelberger et al. 2005a) and fit 
r\ = 1.08 ± 0.17 h^ 1 comoving Mpc. The break 
scale in redshift-space space is larger than the break 



required by the transverse C IV correlation func- 
tion. We attribute the difference to random, pecu- 
liar velocities that impact the line-of-sight correla- 
tion function but have no affect on the transverse 
correlation function. 

• We parameterized the line-of-sight velocity distri- 
bution for the enriched regions and fitted the line- 
of-sight correlation function with ro = 4.0 h^ 1 co- 
moving Mpc ri = 0.42 hr 1 comoving Mpc. The re- 
sulting root-mean-square peculiar velocity was a — 
120 km s _1 with a 3ct upper limit of 300 km s _1 . 
Larger velocities introduced more power at separa- 
tions Av > a than we measure. These speeds may 
entirely reflect the relative motions of galaxies. Be- 
cause the scale r\ is somewhat larger than the virial 
radius, virial motion within galaxies seem less likely 
to dominate. Most importantly, however, we con- 
clude that the typical speeds of outflows over the 
enriched regions are < 200 km s _1 . The maximum 
outflow speed is even lower if the correlation length 
ro is larger than assumed. 

• We measure a higher clustering amplitude for 
C IV absorbers at a median redshift of 2.5 than 
we do at redshift 3.3. This growth in clustering 
amplitude with cosmic time is consistent with the 
gravitational collapse of large structures. The mea- 
surement uncertainties allow on-going enrichment 
however. 

Common C IV systems in spectra of binary quasars re- 
quire enriched regions at least 0.42 h^ 1 comoving Mpc 
in radius. This scale for enriched regions is larger than 
impact parameters at which Steidel etal. (2010) find 
C IV absorption. We attribute the larger size to our 
higher sensitivity to weak C IV systems. Small, but sig- 
nificant, redshift-space distortions explain why the line- 
of-sight correlation function flattens on a larger scale. If 
outflows contribute significantly to these peculiar veloci- 
ties, the outflow velocity could be as high as 200 km s _1 . 
At z ~ 3, our smallest fitted size and largest fitted ve- 
locity give a minimum timescale to disperse the observed 
metals. This timescale exceeds exceeds the median age 
of a Lyman-break galaxy. 

In our qualitative picture, galactic winds disperse met- 
als out to a radius where the outflow stalls or falls back 
towards the galaxy. If these bubbles of metals do not 
overlap, then the size of an enriched region measures the 
turn-around radius of the galactic wind. Large bubbles of 
metals can be made consistent with lower outflow speeds 
by placing the enrichment at a considerably earlier epoch. 
For example, bubbles blown by galactic winds at z> 6, 
need only reach a physical radius of 86 kpc to spread 
metals over 0.42 ft, -1 comoving Mpc by redshift 3. The 
expansion of the universe aids the dispersal of metals in 
this picture. 

The correlation length of the C IV systems only re- 
quires the source of the metals be highly biased like 
Lyman-break galaxies (Scannapieco 2005; Scannapieco 
etal. 2006a). Lower mass galaxies at higher redshift are 
naturally highly clustered (Scannapieco 2005). The lower 
halo mass of early galaxies facilitates the escape of met- 
als (Dekel & Silk 1986; Martin 1999). Their large num- 
ber density may lead to bubble coalescence enlarging the 
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scale of metal-enriched regions well beyond that of an 
individual bubble. 

If the metals detected in these C IV systems were not 
dispersed by galaxies at z ~ 3, then how much earlier 
did dispersal occur? The sizes of the enriched regions 
measured here provide a new reality check for cosmolog- 
ical simulations that address cosmic chemical evolution 
(e.g., Wiersma etal. 2009; Oppenheimer & Dave 2008; 
Kobayashi etal. 2007; Scannapieco etal. 2006b). Some 
simple arguments, however, suggest z> 4.3. Consider 
galaxies at z = 4.3 with the same bias as Lyman-break 
galaxies. The halo mass of such galaxies would be about 
0.2 times that of a typical LBG. The comoving space 
density of such halos is about 5 times higher than that 
of LBGs. In the 723 Myr between redshift 4.3 and 3.0, 
an outflow traveling at 200 km s _1 would reach 0.55 hr 1 
comoving Mpc. Some of these galaxies would still be 
around at redshift 3 and recognizable by their older stel- 
lar population and low mass. (The expansion rate of the 
wind exceeds the cosmic expansion rate, so the wind is 
unlikely to stall.) Taking a cross section of wrf , we esti- 
mate the density of C IV systems (down to our sensitivity 
of 20 mA ) would be dN/dz ~ 5. This value is compara- 
ble to our data if there are 2 to 3 C IV system per galaxy. 
We obtain higher dN/dz, and therefore more systems per 
galaxy if we place enrichment at higher redshift. Only of 
order 10% of of the C IV population is easily attributed 
to enrichment at redshift 3 by Lyman-break galaxies. 

The results desribed above illustrate how quasar to- 
mography can be used to map the distribution of inter- 
galactic metals. The sightlines used in this paper were 
chosen for the relatively small transverse separation of 
the binary quasars. It is possible to observe binaries at 
wider separation. When the correlation of metal- line sys- 
tems is measured over a very broader range in scale, it 
will be possible to simultaneously fit r\ , ro, and 7 to the 
transverse correlation function. The feasibility of mak- 
ing sensitive observations along many pairs of sightlines 
will grow enormously with the next generation of giant 
ground-based telescopes. Our results suggest that mea- 
suring the evolution in the size of enriched regions with 
redshift would directly reveal the IGM enrichment his- 
tory. 
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APPENDIX 

APPENDIX: COMPUTING THE VARIANCE OF OUR CORRELATION MEASUREMENTS 

Single Sightlines 

First we consider measurements along single lines of sight. Numbering all the sightlines in the sample with an index 
£, we can compute the variance in the correlation function for a bin k, representing separations from v — Aw/2 to 
v + Aw/2 as 

a 2 (6) = —^— -xJ2- 2 ((DD)lh (Al) 

Isi { RK )k i 

where (DD) k is the measured number of pairs in sightline £ and bin k, (RR) k is the expected number of pairs in 
sightline £ and bin k for a random distribution, such that ^2 £ (RR) k is the total number of expected pairs in bin k in 
the full data set, and a 2 ((DD) k ) is the variance in {DD) k . 

From Poisson statistics a 2 ((DD) k ) = (DD) k , however as discussed by Mo et al. (1992), this is an underestimate of 
the full variance. To compute the full variance of a single sightline we divide it into many small redshift bins with an 
index i. For simplicity we drop the superscript £ and let Ni be the number of absorbers in one of these bins, where the 
bins are small so that Ni is always zero or 1. To compute the number of pairs in a given sightline, one must compute 
(DD) k = | J2i j NiNjUij where Uij is 1 if the velocity difference between line of sight bins i and j is greater than 

v — Aw/2 and less than v + Aw/2, and Uij = otherwise. 
If we square this we can get an estimate of the variance 

J J2 N i N j N m N n U iij U min , (A2) 

Here there are three possibilities, one is that i, j, in, and n are all different, one is that either i or j falls in the same 
bin as either in or n, and one is that both i and j are matched with m and n. This gives us three terms. If all the bins 
arc different, this gives the average correlation function squared, which does not contribute to the variance. If there 
are two matches between bins, this gives a term that goes as the number of pairs, which leads to the Poisson term. 
Finally, if there is only one match, we obtain a term of the form 

J2 {NiNjN m )U^Ui, m . (A3) 

There is no factor of 1/4 here because there are 4 possible pairings that contribute to this term. The term is equal to 
the joint probability of finding three absorbers: two of them at a distance w ± Aw/2 from the first one. 

There are two possibilities for this distribution. One in which absorption-line systems j and m are on opposite sides 
of absorber i. This possibility contributes 

N ahs x {2nAv[l + £(«)]} x {nAv[l + £(«)]} , (A4) 

where N a bs is the total number of absorption-line systems in a given sightline and n is the average number of absorption- 
line systems per unit length. There is a two on the second term because absorber j can either lie at higher or lower 
redshift than absorber i, but this factor does not appear in the third term because we are assuming that j and m are on 
opposite sides. Note that we are also ignoring the 3-point correlation function, the (usually small) excess probability 
of any absorber being located in bin k above due to both i and j. 

The second possibility for eq. (|A3|) is that absorption-line systems j and in are on the same side of absorber i. In this 
case it is more important that j and m are near each other than at a distance v from absorber i. So this contribution 
looks like 

N ahs x {2nAw[l + £(f)]} x {nA«[l + !(»)]} , (A5) 

where we define 

f(A«) - Av-'' I dv 1 / dv"£(v' - v") , (A6) 



where the double over bar indicates an average over two positions, all positions in which absorption-line systems j and 
m could lie within a bin of size Aw. Again, we have ignored the contribution from the 3 point function, now considering 
only the excess probability due to the fact that j and m are in the same bin. 
Summing up these two contributions and factoring out (DD) k where 

(DD) k = ^{2nAv[l + Z(v)}}, (A7) 



gives 
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C(») + |(Ai;)" 



J2 (NiNjNm) UijUi, m = 2 (DD) k 2nAv 



1 



such that the total variance in sightline £ is given by 



* 2 ((DD)l) = (DD) 



t , 4((DD)l) 2 



N, 



abs 



£(Av)-S(v) 
2[1 + £(«)] 
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(A8) 



(A9) 



In the case in which the correlation function is small, the term in the brackets is just 1, recovering the extra contribution 
to the variance as computed in Mo et al. (1992); but if we know £(i>), then we can also evaluate the term in brackets. 
This variance can be reduced by working with the Landy and Szalay estimator for the correlation function, which 
for a single sightline is given by DD ~ 2 ^^ +RR ■ To compute the variance of this estimator we have to consider two 
additional terms. The first of these arises from the cross-correlation between DD and DR takes the form 

-2 J2 N i N j N m R n U i>j U m , n , (A10) 

i,j,k,m 

where there is a 2 because only the (DD) term is defined with a 1/2 in front. Again we have a term in which i, j, m, 
and n are all different which does not contribute to the variance, and a term in which i and m or j and m are paired, 
which gives 

-4J2 ( N > N 3 R n) UijUi, n . (All) 

i,j,n 

When appropriately normalized to account for the number of absorption-line systems in the random data set this 

- 2A abs [2nAv(l + £(v)]2nAv = -4 (DD) k 2nAv. (A12) 

Finally from the term 4(DRDR) we have a term in which iVj and Nk are pared which, when normalized, gives 

2nAv 



N ahs [2nAvY = 2(DD) t 



! + £(«)' 



(A13) 



Putting this together we get that for the Landy and Szalay estimator, for a single sightline, the total variance is equal 
to the Poisson term, plus a second a term 



2(DD) k 2nAv 



■1 



1 



! + £(«) 



+ 



av)+$(A V ) 



such that the variance in sightline £ is given by 



a'((DD)l) = (DD)l + 



t , 4«DD)«) 2 



JV, 



abs 



m 



l + $(v) 



£(Av) - g v ) 

2(1 + £(«)) 



(A14) 



(A15) 



We do not know the correlation function to evaluate equations (|A9[) and (JA15I) , but as an upper limit on the variance 
we can take the galaxy-galaxy correlation function, taking into account the resolution of our data. This is a power law 
with a slope 7 and correlation length vq = H(z)tq, above w rcs = 50 km/s, and below that. That is 

£(v) = (v/v Q )~\ (A16) 

if v > v rcs and £(v) = otherwise, where rg and 7 are taken to fit the galaxy-galaxy correlation function. In this case 



VO 



if Av > v ICS and £Av = otherwise. 



(2- 7 )(l-7) 



^rcs 
VQ 



^rcs 

Av 



(2-7) Av 1 - 7 



(A17) 



Pairs of Sightlines 



Moving on to the transverse correlation function between two sightlines A and B that make up a pair £, the simplest 
estimator is 

(A18) 



^-ots-tS^' 



where the sum is over all pairs £ whose separations put them in a traverse bin k. The case we are interested in is the 

one in which (D D B ) includes all absorption-line systems that are separated by km/s up to some maximum value, 
Av. 
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Let us again consider the variance for a single sightline. Following our notation from the previous section, we define 
N^ s and N B hs as the number of absorption-line systems in sightline A and B in this pair, and n A and n B as the 
number density of absorption- line systems in sightline A and B. In this case 



(D A D B ) = N A S 2Avn B [1 + fa, A)] = N B hs 2Avn A [1 + fa, At;)], 



(A19) 



where £(r p , Av) = -^ J dv'£(r p ,v'), such that the single bar denotes a single velocity average of £(r p , v 1 ) the corre- 
lation function between two absorption-line systems at a transverse comoving distance of r p and comoving separation 
oft/. 

To compute the variance of a single pair of sightlines, we again divide both sightlines into many small bins and 
square this gives us 

J2 N A N B N A N B U h3 U m , n . (A20) 

This leads to a term where i, j, m and n are all different, which gives the mean value squared, a term where i = m 
and j — n, which corresponds to the Poisson term, and a new term 



J2 (N A N A N, B ) Ui, n U m , n + J2 (N A N B N B ) U hm U, 



j,m, 



(A21) 



Focusing on the first of the two sums in eq. (JA21I) , we can simplify this by considering only the case in which i is at 
a lower redshift than n, and then multiplying by 2. This gives 

53 (N A N A N B ) Ui, m U ntm = N B hs {2n A Av [l + fa, Av)] } [2n A Av [l + fa, Av)] } (A22) 



where 



fa,Av) 



At 



..„..,. dv'[l + fa,v')} dv"{max[fa,v"),m\v"-v'\)}} 

i 1 + ?(?>, Av)\2(Avy Jo i-A„ 



Factoring out (D A D B ) this gives 

J2 (N A N A N B ) Ui, m U n , m = (D A D B ) e 2n A Av [l + fa, Av) 

i,n,m 

Computing the second term in eq. (1A21[) in a similar manner we get 

53 (N A N B N B ) U^U hm = (D A D B ) e 2n B Av [l + fa, Av) 



(A23) 



(A24) 



(A25) 



Thus the variance for pair £ is 



1 1 



a 2 ((D 4 D B )^ = (D A D B Y + ([D A D B Y) 
In the case of the Landy and Szalay estimator, which now goes as D D ~ D R ~ R D +R R 



1 , fa,Av)-fa,Av) 
l + fa,Av) 



r a ru 



(A26) 
, two new cross terms 



arise. The first is a term of the form -2 ({D A D B D A R B ) + (D A D B R A D B )) which leads to 



- 2 53 N A N?RiU id U Jtm + ^ N A R B N B U^U 3 , n = 2 (D A D B ) £ 2Av(n A + n B ). 

\i,j,m i,j,n J 

The second is a term (D A R B D A R B ) + (R A D B R A D B ) , which leads to 

Y^RtRiN B U %] U hm+ 53 N A R B N B U hjUj J = (D A D B ) l 2Av ^ A + n 2 . 

i,j 7 7~n i,j,m,n J s v p? / 

Summing these together then gives 

a 2 ((D A D B )^ = (D A D B ) 1 + ((D A D B )' 



(A27) 



1 1 



fa,Av) \ fa,Av)-fa,Av) 



l + fa,^v) 



1 + fa, Av) 



(A28) 



(A29) 
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In the case of our data, we can estimate £(r p , Av) and £(r p , Av) by carrying out appropriate integrals over 

, n 2 " 

£,(r pi v) 



- 7 /2 
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(A30) 



if v > v rcs or r p =/= and £(r p , v) = otherwise, where II = v/H(z) and ro and 7 are taken to fit the galaxy-galaxy 
correlation function. 
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TABLE 1 
ESI Spectroscopy 



Paired Quasars 


z QSO 


r 


9 


SNR 


Date Observed 


T 




Redshift 


(AB) 


(") 






(s) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SDSSJ0117+3153B 


2.640 


19.86 


11.3 


29.0 


2005 Nov. 


17280 


SDSSJ0117+3153A 


2.625 


20.52 


" 


18.5 


2005 Nov. 


17280 


SDSSJ0225+0048C 


2.700 


20.40 


100.3, 75.3 


12.8 


2006 Nov., 2008 Jan. 


6771 


SDSSJ0225+0048A 


2.820 


20.54 


100.3, 27.4 


12.6 


2005 Nov. 


7100 


SDSSJ0225+0048B 


2.757 


20.40 


75.3, 27.4 


16.5 


2008 Jan. 


7260 


SDSSJ0245-0113A 


2.462 


19.63 


4.5 


20.2 


2006 Nov. 


12400 


SDSSJ0245-0113B 


2.459 


20.51 


" 


16.5 


2006 Nov. 


12400 


SDSSJ0818+0719A 


4.618 


19.66 


132.9 


19.8 


2006 Nov. 


1800 


SDSSJ0818+0718B 


4.178 


20.08 


" 


19.8 


2008 Jan. 


1800 


SDSSJ0925+3859A 


3.135 


19.32 


57.1 


22.5 


2005 Nov. 


6105 


SDSSJ0925+3859B 


2.823 


20.11 


" 


20.8 


2005 Nov. 


10800 


SDSSJ0956+2643A 


3.083 


19.28 


16.5 


19.0 


2007 Apr. 


2400 


SDSSJ0956+2643B 


3.083 


20.52 


" 


14.8 


2007 Apr. 


14740 


SDSSJ0959+1033A 


4.021 


19.34 


44.1 


33.8 


2005 Nov. 


12500 


SDSSJ0959+1033B 


4.021 


20.11 


" 


11.2 


2006 Nov. 


6300 


SDSSJ1021+1112A 


3.830 


20.60 


7.6 


45.7 


2007 Apr. 


21600 


SDSSJ1021+1112B 


3.830 


20.71 


" 


24.8 


2007 Apr. 


21600 


SDSSJ1026+4614A 


3.421 


19.73 


37.1 


24.0 


2008 Jan. 


3600 


SDSSJ1026+4614B 


3.345 


19.73 


" 


17.0 


2008 Jan. 


3600 


SDSSJ1116+4118A 


3.000 


18.18 


13.8 


33.3 


2006 Mar. 


3291 


SDSSJ1116+4118B 


3.000 


19.19 


" 


18.0 


2006 Mar. 


3291 


SDSSJ1144+0959A 


3.146 


18.33 


23.5 


33.6 


2007 Apr. 


2400 


SDSSJ1144+0959B 


2.974 


19.87 




23.1 


2007 Apr. 


10800 


SDSSJ1228+0232A 


3.138 


20.42 


47.3 


11.7 


2008 Jan. 


3600 


SDSSJ1228+0232B 


3.095 


19.83 


" 


11.3 


2008 Jim. 


2400 


SDSSJ1248+1957A 


3.876 


20.23 


64.8 


14.6 


2008 Mar., 2008 June 


4800 


SDSSJ1248+1957B 


3.817 


19.55 


" 


16.5 


2008 Mar., 2008 June 


3000 


SDSSJ1307+0422A 


3.024 


18.05 


8.2 


55.3 


2007 Apr. 


3600 


SDSSJ1307+0422B 


3.024 


19.33 


" 


22.3 


2007 Apr. 


3600 


SDSSJ1314+2818A 


4.821 


21.27 


116.6 


18.6 


2008 July 


2400 


SDSSJ1314+2818B 


4.390 


20.92 


" 


7.5 


2008 July 


900 


SDSSJ1353+4852A 


3.848 


20.43 


37.1 


10.1 


2008 July 


2400 


SDSSJ1353+4852B 


3.848 


20.91 


" 


15.2 


2008 July 


3600 


SDSSJ1420+2830A 


4.312 


20.33 


10.9 


20.8 


2008 Jun. 


10800 


SDSSJ1420+2830B 


4.290 


20.77 


" 


20.5 


2008 Jun. 


10800 


SDSSJ1513-0131A 


3.224 


19.23 


48.1 


22.0 


2007 Apr. 


4200 


SDSSJ1513-0131B 


3.191 


19.91 


" 


18.1 


2007 Apr. 


7200 


SDSSJ1627+4606A 


4.110 


20.66 


34.1 


20.4 


2007 Apr. 


10800 


SDSSJ1627+4606B 


3.812 


20.15 


" 


23.1 


2007 Apr. 


9600 


SDSSJ1054+0216A 


3.976 


19.53 


88.2 


16.2 


2007 Apr. 


2400 


SDSSJ1054+0216B 


3.976 


18.81 


" 


19.8 


2007 Apr. 


2400 


SDSSJ1215-0309A 


4.001 


20.20 


47.1 


13.3 


2007 Apr. 


2400 


SDSSJ1215-0309B 


3.672 


19.81 


" 


15.4 


2007 Apr. 


2400 


SDSSJ1404+4005A 


4.027 


19.62 


47.3 


17.6 


2007 Apr. 


2556 


SDSSJ1404+4005B 


4.027 


21.01 


" 


13.3 


2007 Apr. 


7200 


SDSSJ1541+2702A 


3.621 


20.62 


6.4 


14.1 


2008 Jan., 2008 Mar. 


13782 


SDSSJ1541+2702B 


3.626 


20.71 


" 


10.4 


2008 Jan., 2008 Mar. 


13782 


SDSSJ1542+1733A 


3.261 


20.06 


27.4 


15.2 


2008 June 


3000 


SDSSJ1542+1733B 


2.782 


18.87 


" 


18.5 


2008 June 


1200 


SDSSJ1622+0702A 


3.256 


20.33 


5.8 


14.9 


2008 July 


3000 


SDSSJ1622+0702B 


3.260 


17.19 


" 


86.5 


2008 July 


3000 


SDSSJ1627+2215B 


3.251 


20.55 


35.6 


11.4 


2008 June 


5400 


SDSSJ1627+2215A 


3.709 


19.29 


" 


12.7 


2008 June 


3000 


SDSSJ2157+0015A 


2.560 


20.55 


20.8 


12.7 


2006 Nov. 


4257 


SDSSJ2157+0015B 


2.542 


19.39 


" 


11.6 


2006 Nov. 


3000 





Note. — (1) Quasar name. (2) Quasar redshift. As decribed 
in Hennawi etal. (2009), we measure the wavelegnths of C IV , 
C III, and Si IV and apply the typical shift found between these 
high-ionization lines and the systemic frame, as traced by the 
Mg II doublet (Shen etal. 2007). (3) Quasar r-band magnitude. (4) 
Angular separation of pair. Two values indicate the triplet system. 
(5) Median continuum S/N ratio in the bandpass between Lyct and 
C IV ; the continuum within 5000 km s _1 of cither line is excluded 
from the distribution. (6) Date Observed. (7) Total exposure time. 



